An Enhanced Nonlinear Critical Gradient for Electron Turbulent Transport due to 

Reversed Magnetic Shear 
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The first nonlinear gyrokinetic simulations of electron internal transport barriers (e-ITBs) in 
the National Spherical Torus Experiment show that reversed magnetic shear can suppress ther- 
mal transport by increasing the nonlinear critical gradient for electron-temperature-gradient-driven 
turbulence to three times its linear critical value. An interesting feature of this turbulence is non- 
linearly driven off-midplane radial streamers. This work reinforces the experimental observation 
that magnetic shear is likely an effective way of triggering and sustaining e-ITBs in magnetic fusion 
devices. 
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We aim to investigate, via simulation, recent experi- 
mental evidence on the National Spherical Torus Exper- 
iment (NSTX) [l| of reduced electron thermal transport 
during reversed magnetic shear discharges [2j. 

Understanding and mitigating thermal transport is 
paramount to the development of efficient and econom- 
ical magnetic fusion reactors. In NSTX, turbulent elec- 
tron thermal flux can be very large, often serving as the 
experiment's dominant loss mechanism 0, Yet, by 
reversing the magnetic shear (as defined by s < 0, with 
s = (r/q)dq/dr, r the machine's minor radial coordi- 
nate and q its safety factor), NSTX can reach electron 
temperature gradients that far exceed the limits imposed 
by the onset of unstable electron-temperature-gradicnt- 
driven (ETG) turbulence The present work inves- 
tigates this phenomenon through the first nonlinear gy- 
rokinetic simulations of strongly reversed shear NSTX 
plasmas. We find that reversed magnetic shear, even in 
the absence of E x B shear, suppresses electron turbu- 
lence by greatly enhancing the critical temperature gra- 
dient required for significant thermal transport, and can 
explain observed temperature gradients that significantly 
exceed the limits imposed by linear ETG turbulence. 

Core NSTX plasmas with monotonic q-profiles exhibit 
"stiff profiles": once the electron temperature gradient 
exceeds a critical value, turbulent heat flux relaxes the 
temperature profile back below the critical gradient [f|. 
In such cases, high-k density fluctuation measurements 
show electron-scale turbulence 0- Additional electron 
heating cannot push the temperature gradient greatly 
beyond the linear threshold for ETG drive. Instead, the 
plasma profiles hover near marginal stability, with ETG 
turbulence limiting plasma performance. 

However NSTX can violate the stiff profile phe- 
nomenon. Reversing the magnetic shear triggers an 
"electron internal transport barrier" (e-ITB) in NSTX 
0. Interior to the e-ITB, electron temperatures can grow 



to very high values, while the temperature gradient at the 
e-ITB can reach values two to three times greater than 
the ETG linear critical gradient. Additionally, turbulent 
fluctuations and thermal transport drop [2|, |8[ . An NSTX 
plasma profile with an e-ITB is not stiff around the linear 
ETG critical gradient, but instead at much higher values. 

The present study asks whether a nonlinearly en- 
hanced critical gradient can explain NSTX reversed shear 
c-ITBs. A nonlinear upshift of the critical gradient 
exists for ion-temperature-gradient-driven (ITG) turbu- 
lence 0-11 1, which can allow the normalized ion tem- 
perature gradient to increase by roughly 50% prior to 
the onset of turbulence. While early simulations of ETG 
turbulence showed low transport for s < 12, 13], these 
studies found no enhanced critical gradient for the onset 
of turbulence; turbulent transport remained low for all 
tested gradients. However, these works employed param- 
eters not like those of NSTX. Besides aspect ratio, NSTX 
can have much higher gradients and much lower values of 
8 than those explored in Refs. 12, 13j]. In this context, we 
present the first nonlinear electron-scale gyrokinetic sim- 
ulations of NSTX e-ITB plasmas, using multiple driving 
temperature gradients and values of magnetic shear. 

Herein, we find that reversed magnetic shear, without 
E x B shear, can nonlinearly reduce electron thermal 
transport up to a critical temperature gradient that is 
many times larger than the linear criticalgradient, con- 
sistent with experimental observations @, Q ■ We find the 
effect of reversed shear on ETG turbulence to be much 
more drastic than found that for ITG: the nonlinear crit- 
ical gradient can reach three times the linear threshold, 
consistent with observations of NSTX e-ITBs. 

Our simulations are based upon TRANSP 



14| anal- 



ysis of NSTX shot number 129354 at 232 ms, during a 
radio-frequency-heated e-ITB phase. Experimental de- 
tails are found in Refs. 11 and g. We use GYRO [Tim 
a flux tube geometry j 1 61 ] with a Miller equilibrium 17[ 
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P2c s /a; and the gyroBohm energy flux, 
T e (p s /a) 2 . n e/i , T e/i , m e/D , c and B u 



at the location of the e-ITB. Local parameters include: 
r/a = 0.3, R/a = 1.55, A = -0.27, k = 1.76, s K = 
-0.35, 6 = 0.11, s s = -0.073, C = -2.0 x 10~ 3 , s c = 
0.13, q = 2.40, s = -2.40, p+ = 4.6 x 10" 2 , Z eff = 3.39. 
While Ref. [15j contains full variable definitions, some 
key elements are: the device major radius, R; the full 
minor radius, a; the ion sound speed, c s = ^jT e /mD', 
the ion sound radius, p s — c s {m,DC/ 'eB un it)', the elec- 
tron thermal velocity, v e — \jT e /m e ; the electron gyro- 
radius, p e — v e (m e c/ eB un it); the gyroBohm unit diffu- 
sivity, XGB,i 
QcB,i = n e c t 

represent the background electron/ion density, tempera- 
ture, the electron/deuterium mass, the speed of light and 
the effective magnetic field strength, respectively. Gradi- 
ent lengths are given by the logarithmic derivative with 
respect to the minor radius: Lt = (—T)/(dT/dr). 

The primary ion is cold, low-density deuterium, 
treated gyrokinetically: rii/n e = 0.39, Ti/T e — 
0.22, R/L n% = 1.43, R/L Ti = 2.65 . We use a realistic 
electron mass ratio, p e — y/mi/m e — 60, and the elec- 
tron density scale length remains fixed at R/L He = 1.79. 
The experimental value of the electron temperature gra- 
dient is R/Lt c — 21.5 ± 5. We retain electron-ion col- 
lisions, with v e i — 0.016(a/c s ), but neglect ion-ion col- 
lisions, background E x B flow shear, and finite Debye 
length effects. Our simulations are all electrostatic. 

We follow 24 toroidal modes with 256 radial grid points 
in a box with radial and binormal directions measuring 
L x xL y = 4.26 x 2Ap s = 255.6 x \AAp e , up to a maximum 
kgp e of 1.004, (kgp s ) max — 60.25. Radial convergence 
required Air. < p e . I n ag reement with linear eigenmodc 



Electron Energy Diffusion vs. Temperature Gradient 



simulations of NSTX [18[ , we needed increased resolution 
in other dimensions to reach nonlinear flux con verg ence. 



Details on parameters can be found in Refs. |15l.ll9j. but 
briefly: n T — 44 (the number of mesh points along a 
trapped particle orbit), — 12 (the number of paral- 
lel finite elements), n\ = 12 (the number of pitch angles 
for each sign of parallel velocity, split evenly between 
trapped and passing particles), n t — 12 (the number of 
energy grid points) and e* = 6.0 (the maximum simu- 
lated dimensionless energy). Temporal stability and con- 
vergence required time steps of At (c s /a) = 3.3 x 10~ 4 , 
and simulation times of t max (c s /a) — 10 — 50. Each sim- 
ulation costs 80,000-100,000 CPU hours on the ORNL 
Cray XT system, totaling in all over 5 million CPU hours. 

To probe the existence of a nonlinearly enhanced 
threshold for the onset of electron transport, we scan 
R/Lt b from just above the linear critical gradient, of 
roughly 6, well past physically-relevant values, to over 
25, for three different values of magnetic shear, which 
roughly cover NSTX's reversed shear operating space: 
s = —2.4 (base parameter), —1.0 and —0.2. 

Primarly, we find a nonlinear critical gradient for ETC 
thermal transport that is far larger than the linear crit- 
ical gradient. For the baseline case of s = —2.4, the lin- 




FIG. 1: (color online) Electron heat diffusivities for different 
driving gradients and three values of magnetic shear. The 
horizontal line at Xe = 1 m 2 /s indicates experimentally rele- 
vant levels. 



Shear Dependence of Critical Gradients 
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FIG. 2: (color online) Critical gradients as a function of mag- 
netic shear. ETG becomes linearly unstable at gradients 
above z c (solid with square). Above z^ L (dashed with cir- 
cles), turbulent Xe exceeds experimental values, 1 m 2 /s. The 
dotted line with diamonds represents the nonlinear upshift of 
the critical gradient, Az = z^ L — z c . For comparison, the 
original cyclone ITG test case found an upshift in the critical 
gradient for transport that extended R/L^ by 2, from z c — 4 



to Zc L = 6 at s 
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ear critical gradient for turbulence is where the growth 
rate becomes positive, z c = (R/ Lr e ) critical, linear ~ 6. 
Yet the heat flux remains low for driving gradients be- 
low z»l = (R/L T J crMcal ^ nonlmear a 19, the effective 
nonlinear critical gradient for transport. The difference, 



Az 



z c , represents the nonlinear upshift of the 



critical gradient for transport. In this case, Az « 13. 

The strength of this upshift depends upon magnetic 
shear. Figure[T]shows electron heat diffusivities as a func- 
tion of R/Lx e and s. The horizontal dashed line corre- 
sponds to experimentally relevant levels, at % e = 1 m 2 /s. 
We define z^ L as the gradient above which \& exceeds 
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h e s = -2.4 R/L Te = 9.3 



Time-Averaged Heat Flux vs. Toroidal Mode 
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FIG. 3: (color online) Poloidal cross sections of saturated 
density fluctuations below and above z^ L for s = —2.4. 



this value, and plot this gradient and the corresponding 
values of z c for the cases in Fig. [5] As s decreases, z c , 
and Az increase, up to very large values: for the 



baseline value of s — — 2.4, z 
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The nature of the saturated plasma turbulence changes 
drastically above and below z~? L . Figure [3] shows a 
poloidal cross section of computed electron density fluc- 
tuations above and below z^ L . For clarity we show 
only |n e | < 2%. The top of Fig. [3] shows that below 
z^ L , fluctuations appear on the midplane. Their ampli- 
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FIG. 4: (color online) Time-averaged heat flux spectra for two 
values of (kep a ) ma x, at R/Lt c = 21.8, and the fastest growing 
linear mode. t a v g = [4.0,8.0] (a/c s ). 



Modal and Poloidal Dependence of Density Fluctuations 




FIG. 5: Density fluctuations as a function of 6 and kgp a 
showing off- midplane peaking in kep s = 13. 



tude is small {{n e ) rms » 0.3%, \n e \ max ~ 1.3%), and 
the sheared structure of the streamers is clearly visi- 
ble. Eddies quickly rotate towards the vertical, thereby 
limiting their effect on transport, consistent with ear- 
lier investigations of reversed magnetic shear [3, l20l |. 
Above z^ L , a very different picture emerges, as in the 
top of Fig. |3l Not only are these fluctuations stronger 
((n e ) rms ~ 1-1%, \n e \ max w 6.0%), their behavior is 
drastically different. Streamers appear out of the top and 
bottom of the flux tube, instead of along the midplane. 

The off-midplane streamers are localized to a nar- 
row spectral band and appear to be nonlinearly driven 
by the peak of the linear growth spectrum. Figure |4] 
shows the linear spectrum and time-averaged nonlinear 
heat flux per toroidal mode for two different values of 
(kgp s )max- Though the fastest growing mode is cen- 
tered near kgp s = 40, the late-time heat flux consists 
of a broad spectrum about kgp s = 23 and a strong peak, 
near kgp s — 13. The poloidal dependence of density fluc- 
tuations in Fig. [5] shows that, while the broad spectrum 
centers on the midplane, the spike is the heat flux from 
the off-midplane streamers. By reducing k max , this heat 
flux drops, while the background spectrum is modestly 
affected. This implies not only nonlinear energy cascade 
to the broad background turbulence, but also preferen- 
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tial coupling to the off-midplane streamers. We note that 
both the poloidal angle of the off-midplane streamers and 
their position in k-space vary with temperature gradi- 
ent and shear. Above the upshifted nonlinear critical 
gradient, we find that nearly one quarter of the turbu- 
lent transport is caused by a single slowly-growing mode 
driven to high amplitudes by the spectrum of faster grow- 
ing, higher wavenumber modes. 

The enhanced upshift of the critical gradient for ETG 
transport might be explained by secondary instability 
theory 12 , [ljj HH S3 • Primary instabilities grow in the 
plasma until they are balanced by a growth in a sec- 
ondary mode, which serves as an energy sink for the pri- 
mary mode. The growth of the secondary mode depends 
not only upon the properties of the primary mode (such 
as amplitude and wavenumber), but also upon some 
properties of the background plasma, such as R/Lt c and 
s. Gyrofluid estimates [ij} suggest negative magnetic 
shear reduces ETG transport by enhancing the damping 
due to the secondary mode. While this motivates the 
reduction in transport in these NSTX cases, to fully un- 
derstand this strong upshift in the critical gradient, it 
would be useful to extend the work on gyrokinetic sec- 
ondary theory [22} to reversed shear cases. An interesting 
note is that we, likeRefs. [HEi, 

find low transport for 
adiabatic ion (ETG-ai) simulations. This is the opposite 
behavior than the ETG-ai transport runaway found at 
high positive shear (23|. However, kinetic ions can affect 
the growth rate of the gyrokinetic ETG secondary insta- 
bility [13] • The discovery of an ETG-ki nonlinear critical 
gradient further motivates re-examining this theory in 
the context of reversed shear ETG turbulence. 

In summary, our goal was to numerically investigate 
the presence of e-ITBs during reversed shear NSTX dis- 
charges. Our first-of-a-kind nonlinear ETG simulations 
of NSTX e-ITBs show that reversed magnetic shear, even 
in the absence of background E x B flow shear, is very ef- 
fective at suppressing electron transport by significantly 
enhancing the nonlinear critical gradient for ETG turbu- 
lent flux. The upshifted gradient, which can approach 
three times the linear critical gradient, is much larger 
than predicted upshifts for ITG turbulence 

Above this upshifted threshold, the turbulence takes 
on a different form than that first associated with 
ETG turbulence 12|, ll3j . Characterized instead by off- 
midplane streamers, this variation in structure could 
have implications for high-k density scattering measure- 
ments during these discharges. Interpreting these results 
with a synthetic high-k diagnostic (such as in Ref. (24[) 
and experimentally looking for these structures (by mov- 
ing the high-k measurement location, for instance) , could 
help elucidate the "bursty" turbulence measured in some 
reversed shear NSTX e-ITBs 0. 

These results help explain not only NSTX e-ITBs, but 
also observations on other magnetic fusion devices. The 



RFX-Mod Reversed Field Pinch has found e-ITBs with 
non-monotonic q profile s 1251 . a feature also linked to ITB 



formation in tokamaks 2(3-30]. As stellarators' external 
coils dictate the magnetic topology, it may be possible to 
optimize stellarator designs for e-ITBs. In short, shear 
reversal, even without background flow shear, could allow 
magnetic fusion devices to sustain electron temperature 
gradients much larger than the limits imposed by the 
linear onset of instability. 
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